236 


JOURNAL OF THE ATMOSPHERIC SCIENCES 


Volume 69 


The Impact of Dry Midlevel Air on Hurricane Intensity in Idealized 
Simulations with No Mean Flow 

Scott A. Braun and Jason A. Sippel 

Laboratory for Atmospheres, NASA Goddard Space Flight Center, Greenbelt, Maryland 

David S. Nolan 

Rosenstiel School of Marine and Atmospheric Science, University of Miami, Miami, Florida 


(Manuscript received 24 November 2010, in final form 22 June 2011) 

ABSTRACT 

This study examines the potential negative influences of dry midlevel air on the development of tropical 
cyclones (specifically, its role in enhancing cold downdraft activity and suppressing storm development). The 
Weather Research and Forecasting model is used to construct two sets of idealized simulations of hurricane 
development in environments with different configurations of dry air. The first set of simulations begins with 
dry air located north of the vortex center by distances ranging from 0 to 270 km, whereas the second set of 
simulations begins with dry air completely surrounding the vortex, but with moist envelopes in the vortex core 
ranging in size from 0 to 150 km in radius. 

No impact of the dry air is seen for dry layers located more than 270 km north of the initial vortex center 
(~3 times the initial radius of maximum wind). When the dry air is initially closer to the vortex center, it 
suppresses convective development where it entrains into the storm circulation, leading to increasingly 
asymmetric convection and slower storm development. The presence of dry air throughout the domain, in- 
cluding the vortex center, substantially slows storm development. However, the presence of a moist envelope 
around the vortex center eliminates the deleterious impact on storm intensity. Instead, storm size is signifi- 
cantly reduced. The simulations suggest that dry air slows intensification only when it is located very close to 
the vortex core at early times. When it does slow storm development, it does so primarily by inducing out- 
ward-moving convective asymmetries that temporarily shift latent heating radially outward away from the 
high-vorticity inner core. 


1. Introduction 

High relative humidity (RH) in the middle tropo- 
sphere has long been recognized as an important factor 
in determining where tropical cyclones form (Gray 1975, 
1979, 1998; McBride 1981). Its favorable role was viewed 
more in terms of being a necessary climatological con- 
dition rather than being a determining factor in whether 
or not individual cloud clusters went on to develop into 
tropical cyclones (McBride 1981; McBride and Zehr 
1981). However, DeMaria et al. (2001) showed that their 
formulation of a genesis parameter, of which midlevel 
moisture is a part, can provide some useful information 
as to the probability of tropical storm formation. Kaplan 
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and DeMaria (2003) showed that high values of 850- 
700-hPa relative humidity generally favor rapid inten- 
sification of tropical cyclones. 

Kimball (2006) examined the impact of dry intrusions 
by perturbing initial moisture in simulations of Hurricane 
Danny (1997). Kimball varied both the magnitude of the 
inner-core moisture anomaly (4 g kg -1 variations in peak 
magnitude, maximum in the boundary layer and de- 
creasing with height) and its size (from 250 to 600 km). 
As might be expected, an initial vortex with higher 
moisture content (for fixed size) generally led to more 
intense storms, while more extensive moisture anom- 
alies typically led to increased areal extent of rainbands 
and a larger area of storm-strength winds (17 m s -1 ). 
Kimball claimed that dry air intrusions into systems 
with smaller moist envelopes contributed to weakening 
of those cases, although the differences in minimum 
central sea level pressure were generally less than 5 hPa 
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prior to landfall for the experiments with peak moisture 
of 19 g kg -1 (her A19 simulations). Because moisture 
perturbations extended into the boundary layer, the role 
of midlevel moisture anomalies alone is somewhat am- 
biguous in these simulations. 

Hill and Lackmann (2009) performed idealized sim- 
ulations to examine the impact of environmental moisture 
on storm size. Using a fixed initial vortex and inner-core 
thermodynamic conditions (80% RH within a 100-km 
radius), they varied the environmental relative humidity 
between 20% and 80%. In results that were comparable 
to Kimball (2006), they found that higher environmental 
humidities led to increased outer rainband production, 
larger storms, and broader storm-force wind distribu- 
tions. In terms of storm intensity, differences among the 
20%, 40%, and 60% relative humidity cases were mini- 
mal while the 80% humidity case had a lower minimum 
central pressure but nearly identical maximum winds. 
Their results suggest that the environmental humidity has 
a critical impact on storm size, but a much smaller impact 
on storm peak intensity (as measured by conventional 
parameters such as minimum pressure or maximum wind 
speed). 

Dry midlevel air is also one of the potential inhibiting 
influences of the Saharan air layer (SAL) on the de- 
velopment of tropical cyclones (Dunion and Velden 2004; 
Wu et al. 2006; Jones et al. 2007; Wu 2007; Dunion and 
Marron 2008; Sun et al. 2008, 2009; Shu and Wu 2009). 
Dunion and Velden (2004) suggested that dry SAL air 
can negatively impact tropical cyclones by fostering en- 
hanced cold downdrafts (Emanuel 1989; Powell 1990) 
and lowering the convective available potential energy 
(CAPE) within tropical cyclones. Shu and Wu (2009) 
suggested that storms often weaken when dry air en- 
croaches to within 360 km of the storm center. However, 
in a composite of National Centers for Environmental 
Prediction (NCEP) Global Data Assimilation System 
(GDAS) final (FNL) analyses for 41 storms, Braun (2010) 
found statistically insignificant differences in relative hu- 
midity between subsets of storms that strengthened and 
that weakened in the days after formation, suggesting that 
dry SAL air only weakly impacts intensity once distur- 
bances reach tropical storm strength. In a statistical anal- 
ysis of a 30-member ensemble of simulations of Tropical 
Storm Debby (2006), Sippel et al. (2011) found sensitivity 
of storm intensification to dry SAL air only up to the early 
depression stage of the storm, consistent with the results of 
Braun (2010). 

Many studies have focused on dry air as a key mech- 
anism for hurricane suppression or weakening, as a re- 
sult either of the SAL (Dunion and Velden 2004; Wu 
et al. 2006; Jones et al. 2007; Wu 2007; Shu and Wu 2009; 
Sun et al. 2008, 2009) or of other dry air sources (Kimball 


2006). To investigate this process more carefully, this 
study uses a set of idealized simulations to examine the 
impact of a dry midlevel air layer. These simulations 
also use an environment with no mean flow. More com- 
plex simulations with mean flow, including vertical wind 
shear, are reserved for a future study. Section 2 describes 
the model setup and experiments, section 3 describes 
impacts on storm intensity, section 4 describes impacts on 
storm structure, and section 5 explains how the variations 
in structure relate to the differing intensification rates of 
the storms. Conclusions are given in section 6. 

2. Model setup and experiments 

This study employs the Advanced Research version of 
the Weather Research and Forecasting (WRF) modeling 
system (version 3.1; Skamarock et al. 2005) to conduct 
idealized simulations of the interaction of developing 
tropical cyclones with dry midlevel air. The model has 
been applied successfully to the simulation of real tropi- 
cal cyclones (Davis et al. 2008a,b; Nolan et al. 2009; Vizy 
and Cook 2009; Braun et al. 2010; Wang et al. 2010, 
among many others) and for idealized simulations (Nolan 
2007; Nolan et al. 2007b; Nolan and Rappin 2008; Hill 
and Lackmann 2009). In particular, Davis et al. (2008b) 
showed that the WRF model performance in real-time 
simulations was competitive with, and in some cases 
superior to, other operational models. Here, three grids 
nesting down to 2-km horizontal grid spacing are em- 
ployed in order to reasonably represent the convection. 
The outer grid has a horizontal grid spacing of 18 km 
and contains 240 X 240 grid points in the x and y di- 
rections. Two nested meshes are used with the following 
grid spacings and grid dimensions: 6 km and 120 X 120, 
and 2 km and 240 X 240. All grids use 49 vertical levels 
with a model top at 20 km. Physics options include the 
Yonsei University boundary layer scheme (Noh et al. 
2003; Hong et al. 2006), the fifth-generation Pennsylvania 
State University-National Center for Atmospheric 
Research (NCAR) Mesoscale Model (MM5) similarity 
theory surface-layer scheme (Zhang and Anthes 1982; 
Skamarock et al. 2005), the Kain-Fritsch cumulus scheme 
(Kain and Fritsch 1990, 1993; Skamarock et al. 2005) on 
the 18-km grid, and the WRF single -moment three-class 
simple ice cloud microphysics scheme (Hong et al. 2004) 
on all grids. Radiative processes are not included. 

Initial and boundary conditions are derived loosely 
following Nolan et al. (2007a) and Nolan and Rappin 
(2008). The boundary conditions are doubly periodic. 
The domain is on an / plane with background Coriolis 
parameter of / = 5.0 X 10 -5 s -1 . The initial vortex is 
specified as a modified Rankine vortex with solid body 
rotation out to a radius of maximum winds (RMW) of 
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Table 1. Simulation names and descriptions. 


Simulation 

name 

Description 

CNTL 

Control run with uniform non-SAL 
environment 

DRY270 

Dry air located northward of 270 km north of 
vortex center 

DRY144 

Dry air located northward of 144 km north of 
vortex center 

DRY90 

Dry air located northward of 90 km north of 
vortex center 

DRY0 

Dry air located northward of vortex center 

DRYALL0 

Dry air throughout domain, including vortex 

DRYALL75 

Dry air throughout domain except within 
75-km radius 

DRYALL150 

Dry air throughout domain except within 
150-km radius 


100 km and a wind decay of l/r" beyond the RMW, with 
a = 0.5. The maximum wind speed is 15 m s _1 at 4.5-km 
altitude. The sea surface temperature is set to 28°C 
(301.15 K). The background environment representing 
non-SAL tropical air is specified by the Dunion and Marron 
(2008) non-SAL sounding. For the control run (CNTL) that 
excludes a dry midlevel air layer, this environment is 
specified for the entire domain. In these simulations, 
there is no mean environmental flow, which simplifies 
the problem to the issue of whether the interaction be- 
tween the dry air and the storm circulation is sufficient 
to affect storm evolution. 

Several different configurations for the dry air layer 
are utilized. The first, labeled DRY270, places a dry air 
layer at all grid points farther north than 270 km (15 grid 
points on the outer domain) north of the initial storm 
center. North of this boundary, the relative humidity in 
the 850-600-hPa layer is set to 25%. Note that composite 
relative humidity fields for weakening storms in Fig. 13 of 
Braun (2010) show the dry air typically farther than 
400 km from the storm center, so the position of the dry 
air boundary is closer to the vortex than is suggested by 
NCEP final analyses. Additional simulations, labeled 
DRY144, DRY90, and DRY0 move the dry air 
boundary successively closer to the vortex center until in 
the last case the dry air boundary is at the center. The final 
set of simulations involves the placement of the dry air 
layer over the entire domain. In DRYALL0, the dry air 
extends throughout the vortex and its surroundings. 
Given that storms typically form within an envelope of 
moist air, in DRYALL150 and DRYALL75 a moist 
envelope (using the non-SAL sounding) is prescribed 
from the center of the vortex to a radius of 150 and 75 km, 
respectively, but with dry air otherwise surrounding the 
vortex. A summary of all experiments is given in Table 
1. In the discussion that follows, the minimum sea level 


pressure and maximum instantaneous wind speed are 
taken as measures of storm intensity. 

The storm center was determined, as in Braun (2002), 
at every model output time using the surface pressure 
field. Rather than using the location of the minimum 
pressure, we used the horizontal distribution of pressure 
to determine an approximate geometric centroid. The 
location of the minimum pressure was used as a first 
guess for the center, and a variational approach then 
adjusted the center location to minimize the azimuthal 
variance of the pressure field at all radii within 150 km. 

3. Impacts on storm intensity 

The results from Fig. 1 suggest that the dry midlevel 
air can act as a brake on development, but only if it is 
able to penetrate to very small radius during the early 
stages of development. First, the CNTL run with no dry 
layer starts with an initial minimum central pressure of 
1012 hPa and takes about 2 days to reach 1000 hPa and 
begin a more rapid intensification process. By day 7, the 
central pressure drops to a minimum of 940 hPa and the 
maximum wind speed reaches —50 m s , remaining ap- 
proximately steady thereafter. Meanwhile, the minimum 
pressure and maximum wind speed in DRY270 are almost 
identical to those in the CNTL run, suggesting little impact 
of the dry air on storm intensity for initial separation 
distances greater than 270 km. As the separation distance 
is reduced from 270 km to 144 (DRY144), 90 (DRY90), 
and 0 km (DRY0), the dry air increasingly slows or de- 
lays the intensification of the vortex, although all cases 
reach approximately the same maximum intensity by the 
end of the simulations. 1 

To illustrate why the dry air has little impact on storm 
intensity for the DRY270 case despite its proximity to 
the vortex. Fig. 2 shows the azimuthally averaged rela- 
tive humidity and radial winds for the DRY270 simula- 
tion at 5 days. Dry air with relative humidity less than 
50% does not penetrate (in the average) inside of 300 km 
from the vortex center. The primary reason is that the 
radial flow in the 850-600-hPa layer is characterized by 
weak inflow (<2 m s -1 ) outside —225 km and weak 
outflow at smaller radius. This is similar to several pre- 
vious simulation studies (Liu et al. 1997; Braun 2006; 
Smith et al. 2009) where the bulk of the radial inflow is in 
the surface boundary layer, typically below —900 hPa or 


1 The simulations achieve the same final intensity because, given 
enough time, the dry air is eventually moistened so that all systems 
reach their maximum potential intensity (within the context of the 
modeling system and the choices for parameterizations). The 
maximum potential intensity here is determined by the imposed 
SST and mean sounding. 
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Fig. 1. Time series of (a) minimum sea level pressure and 
(b) maximum wind speed for all simulations. The solid line is the 
CNTL run, the dashed lines are the runs with dry air north of 
a specified latitude, and the dotted lines are for runs with dry air 
throughout the domain except for a moist envelope within the vortex. 


1-km altitude. This figure shows that, in the absence of 
mean environmental flow, there is no mechanism for the 
dry air to penetrate into the inner core of convection and 
that the only way for it to affect storm development is to 
begin in the inner-core region at initial times. 

Figure 3 illustrates the situation where dry air exists in 
the initial inner core. In this figure, relative humidity at 
3-km altitude from the 18-km grid of DRYO is shown for 
selected times. The southern half of the vortex is initially 
collocated with an area of enhanced humidity (Fig. 3a) 
resulting from the cool anomaly below the balanced 
midlevel vortex, and the boundary of the dry midlevel air 
goes through the center of the vortex. The cyclonic winds 
associated with the vortex extend well into the dry air 
and, after 2 days (Fig. 3b), they have wrapped this dry air 
nearly all the way around the vortex. In addition to the 
southward displacement of dry air, the moister air ini- 
tially on the southern side of the vortex gets displaced 
northward into the dry air mass. Over time (Figs. 3c, d), 
the dry air gets increasingly wrapped around the vortex 
and axisymmetrized, lowering the overall relative hu- 
midity in the environment of the vortex compared to 
CNTL. The effect on storm intensity (Fig. 1) is to delay or 
slow intensification of the storm by up to 2 days, although 
the storm reaches a maximum intensity similar to that of 
CNTL by the end of the simulation period (i.e., 8 days). 


As might be expected, the largest impact of dry air is 
found when the dry layer extends across the entire do- 
main (case DRYALLO). In this case, deep convection 
does not begin until the end of the third day, after shallow 
convection has gradually moistened the midlevels, with 
intensification of the vortex then ensuing (Fig. 1). The 
storm begins a period of more rapid intensification on the 
fifth day, approaching the maximum intensity of the other 
simulations after 8 days. 

Both DRYALLO and DRYO are generally unrealistic 
for development of actual tropical disturbances since these 
disturbances virtually always form within some moist en- 
velope of air associated with easterly waves (Dunkerton 
et al. 2009; Hopsch et al. 2010). This fact suggests that 
a more realistic initial condition for an environment with 
dry air surrounding the vortex would include a moist 
envelope with at least moderate (>60%) relative hu- 
midity within some radius from the initial storm center. 
To examine the impact of a moist envelope collocated 
with the vortex, two experiments were performed in 
which the dry midlevel air (from case DRYALLO) within 
some radius R was replaced by the non-SAL sounding. 
Given that the initial radius of maximum wind is 100 km, 
we tested R = 150 km and R = 75 km. Results from the 
latter case will be emphasized since both simulations 
produce similar results (Fig. 1). 

Initially, moist (>80% RH) conditions exist only 
within the small predefined region near the storm 
center and the circulation is otherwise surrounded by 
very dry air. Over time, the moist region expands as 
convection increases winds within the boundary layer, 
driving larger fluxes of sensible and latent heat. The 
intensity of the storm in DRYALL75 as a function of 
time (Fig. 1) is essentially identical to that in the CNTL 
case, suggesting that a vortex with even a modest-sized 
moist envelope will not necessarily be adversely af- 
fected by dry midlevel air, even when completely sur- 
rounded by it. This result is qualitatively consistent 
with those of Hill and Lackmann (2009), who found that 
relative humidity outside a 100-km radius had very lim- 
ited impact on storm central pressure and maximum 
wind speed. 

4. Impacts on storm structure 

The intrusion of dry midlevel air can produce marked 
changes in storm structure. To illustrate the character- 
istics of these changes. Figs. 4 and 5 show horizontal 
distributions of simulated midlevel (3 km) relative hu- 
midity and near-surface radar reflectivity, respectively, 
for 12 h and 3 days into the simulations from CNTL, 
DRYO, and DRYALL75. At 12 h into the simulation, 
the CNTL run shows two nearly concentric rings of high 
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Fig. 2. Azimuthally averaged RH (shading; contours at 30%, 50%, 70%, and 90%) and radial 
velocity (contours, 2 m s _1 intervals, negative values dashed) from the 18-km grid of simulation 
DRY270 at 5 days. 


relative humidity and reflectivity (Figs. 4a and 5a). An 
inner ring lies just inside the initial radius of maximum 
wind and an outer ring is along the edge of a circular cold 
pool region formed by the initial convection (to be dis- 
cussed in section 5). These rings should not be confused 
with convective organization that would be associated 
with secondary eyewall formation; rather, they are a 
consequence of the strong symmetry of the initial con- 
ditions. In contrast, in the DRY0 simulation, dry air has 
wrapped into the eastern portion of its vortex (Fig. 4c), 
largely suppressing any convection there (Fig. 5c). As 
a result, the precipitation is highly asymmetric in the 
outer convective ring and very limited in the inner ring. 
The horizontal precipitation structure of the storm in the 
DRYALL75 simulation (Fig. 5e) shows suppressed de- 
velopment of the outer convective ring but several con- 
vective cells within an inner convective ring. 

The suppressing influence of the dry air is illustrated in 
Fig. 6a, which shows the horizontal distribution of relative 
humidity and simulated radar reflectivity 6 h into the 
DRY0 simulation. The deep convection, as indicated by the 
higher reflectivities, is contained within the moist portion of 
the vortex. Dry air is wrapping around the southern side of 
the vortex. Convection within this dry region is generally 
limited to a depth of just 3 km (Fig. 6b, between —170 and 
300 km) except along the western edge of the cold pool 
(Fig. 6c, near 130-140 km) where it penetrates to just above 
4 km. The curvilinear pocket of moderate to high humidity 
within the broader dry air mass is associated with this 


— 4-km-deep convection along the cold-pool edge. Vertical 
cross sections through this feature (Figs. 6b,c) show that the 
convection generally does not penetrate above the dry 
layer, while convection in the moist layer easily reaches 
10 km or higher. The shallow convection in the dry air re- 
sults from a reduction of boundary layer equivalent po- 
tential temperature 9 e caused by cold downdrafts, thereby 
reducing the parcel 6 e entering cloud base, and the rapid 
entrainment of dry air as the convective plumes rise. 

Generally similar differences between CNTL and 
DRY0 continue over the next several days of the simula- 
tions. By the end of day 3, the CNTL simulation shows 
a very symmetric system of tropical storm strength (Figs. 
4b and 5b), with a nascent eyewall as well as inner and 
outer rainband structures. In the DRY0 simulation, the 
only remaining very dry air (<50%) is found at a radius of 
about 200 km (Fig. 4d). Although a weak asymmetry in 
the relative humidity field exists, with humidities between 
70% and 80% prevalent on the eastern side of the storm 
compared to higher humidities elsewhere, the dry air has 
largely been axisymmetrized. Despite the axisymmetriza- 
tion of the humidity field and formation of a nascent 
eyewall (Fig. 5d), the convection remains highly asym- 
metric, with the major portion of the precipitation in the 
southwestern quadrant. By 3 days into DRYALL75, 
a symmetric eyewall has formed (Fig. 5f), but with a radius 
that is only half that in the CNTL simulation. 

The evolution of the precipitation structure in the 
CNTL, DRY270, DRY144, DRY90, DRY0, and 
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Fig. 3. RH at 3-km altitude from the DRYO simulation (18-km coarse grid domain) at (a) the initial time, and after 
(b) 2, (c) 4, and (d) 6 days. The color scale is shown in (a). 


DRYALL75 simulations is depicted in Fig. 7 in terms 
of the radial and azimuthal variations of the 3-km-level 
rain mixing ratio with time for the first 6 days of the 
simulations. In CNTL (Fig. 7a), convection begins near 
a radius of 70 km and propagates outward along an 
outflow boundary. Two more cycles of convection and 
outflow occur before the precipitation becomes more 
continuous in time after 30 h. The precipitation in the 
nascent eyewall contracts and intensifies with time up 
to 4.5 days and then moves slowly outward between 5 
and 6 days. The precipitation field is very symmetric at 
early stages (Fig. 7b), but it exhibits a weak wavenunrber-1 
asymmetry between 3 and 5 days. 

When dry air is placed 270 km north of the vortex 
center, the evolution of the radial distribution of rain 
(Fig. 7c) and the storm intensity are very similar to 
CNTL. A stronger wavenumber-1 asymmetry (Fig. 7d) 
becomes noticeable as early as 12 h and better defined 


after 30 h. Although the rain field is asymmetric, pre- 
cipitation still generally encloses most of the developing 
eye, even if only weakly. When the dry air is initially 
placed at 144 km from the center, the rain field (Fig. 7e) 
still steadily contracts, but it intensifies more slowly. A 
more pronounced change occurs in the azimuthal vari- 
ation of the rain field (Fig. 7f), which has a much more 
well-defined wavenumber-1 asymmetry. The asymmetry 
is sufficiently strong that a significant portion of the 
forming eyewall has little or no rain through much of 
the 6 days shown. 

When dry air is initially placed at a radius of 90 km, 
just inside the initial 100-km radius of maximum winds, 
an important change in precipitation evolution occurs 
(Fig. 7g). Unlike the fairly continuous contraction of the 
precipitation fields in the earlier cases, the precipitation 
field shows little or no contraction between 1.5 and 
3.5 days. After that time, the rain maximum contracts 
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Fig. 4. RH at 3-km altitude from the 2-km resolution domain for (left) 12 h and (right) 3 days into the simulation, 
for (top) the CNTL simulation, (middle) the DRY0 simulation, and (bottom) the DRYALL75 simulation. Thick 
solid contours are drawn at 50% RH. 
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Fig. 5. Simulated radar reflectivity at 0.5-km altitude from the 2-km resolution domain for (left) 12 h and 
(right) 3 days into the simulation, for (top) the CNTL simulation, (middle) the DRY0 simulation, and (bottom) 
the DRYALL75 simulation. Thick solid contours show the 50% RH contours from Fig. 4. 
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Fig. 6. (a) RH (shading) and simulated radar reflectivity (contours, 10-dBZ interval, starting at 5 dBZ) at 3-km 
altitude from the innermost grid of the DRYO simulation at 6 h. The white line indicates the location of the vertical 
cross sections. (b),(c) Vertical cross sections of (b) RH (20% contour interval) and (c) 8 e (3-K contour interval). The 
color scale for (a) and (b) is in the lower-left part of the figure. The color scale for (c) is in the lower-right portion of 
the figure. 


from about 60-km radius to about 25-km radius during 
the next 24 h and then intensifies after 4.5 days. As 
in the DRY144 simulation, the structure in DRY90 is 
characterized by a pronounced wavenumber-1 asym- 
metry through the first 6 days of the storm’s evolution 
(Fig. 7h), with little or no rain over half the eyewall. 
A very similar evolution is seen in DRYO (Figs. 7i,j), 
with the main difference being a greater reduction 
in rain during the first 2 days of simulation as the dry 
air occupies half the area of the initial vortex. Similar 
to DRY90, precipitation remains weak and steady 
prior to the time of intensification (near 4.75 days), at 
which time the radius of rainfall contracts and the rain 
amount gradually increases. The pronounced wave- 
number-1 asymmetry (see Fig. 5d) is also evident in 
DRYO. 

When a small moist envelope is included, but the 
storm is otherwise surrounded by dry air (DRYALL75; 
Figs. 7k,l), the evolution is very similar to CNTL except 
that rain is concentrated near a radius of 20-30 km 
rather than 30-40 km in CNTL. The radial structure is 
characterized by rapid contraction of the nascent eyewall 
during the first 2 days and slow expansion after 5 days. 
Outward-moving rainbands occur frequently through 
the first 5 days of simulation. The azimuthal structure 
shows a somewhat more pronounced wavenumber-1 
asymmetry compared to CNTL, but this asymmetry is 
much weaker than that found in the other simulations 
with dry air. 


5. The relative impacts of cold pools and enhanced 
storm asymmetry 

The intrusion of dry air into the vortex can have mul- 
tiple effects including 1) suppression of deep convective 
development because of boundary layer cooling from 
downdrafts, 2) suppression of deep convection by en- 
trainment of dry air into rising clouds, and 3) generation 
of precipitation asymmetries that persist well beyond the 
time of axisymmetrization of the humidity field. In this 
section, we examine the roles of these mechanisms in the 
impact of the dry midlevel air on storm intensity. 

Figure 8 shows plots of the near-surface (90 m) d e and 
simulated radar reflectivity at 12 h for the CNTL, DRY90, 
DRYO, and DRYALL75 simulations. The azimuthally 
averaged near-surface 6 e distributions as a function of 
radius and time are shown in Fig. 9 for each of these cases. 
In CNTL (Fig. 8a), a symmetric cold outflow is generated 
by the initial convection, but the interior region where 
surface winds are largest recovers most quickly. A nearly 
identical pattern is found in the DRY270 simulation (not 
shown). Figure 9a shows the formation and outward 
movement of the initial cold pool, as well as several 
additional outward-moving cold pools that occur over 
the next 2 days. After about 2.5 days, d e rapidly in- 
creases in the eye and increases more slowly outside of 
the eyewall as the storm intensifies. 

When dry air is placed 90 km north of the vortex 
center (DRY90), it wraps around the western side of the 
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Fig. 7. (column 1), (column 3) Radius vs time distributions of azimuthally averaged 3-km rainwater mixing ratio. The corresponding 
color scale is the lower-left scale, (column 2), (column 4) Azimuth vs time distributions of the radially averaged (0-150-km radius) 
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panel. 
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Fig. 8. Equivalent potential temperature d e at the lowest model level (90 m) at 12 h into the simulations from the 
2-km grid for (a) the CNTL run, (b) the DRY90 run, (c) the DRY0 run, and (d) the DRYALL75 run. Thick blue 
contours show the 50% RH contours at 3-km altitude. Thin white contours show simulated radar reflectivity at 2-km 
altitude every 10 dBZ starting at 25 dBZ. The thick red X marks the center location. 


vortex and enters the precipitation area near its south- 
ern end by 6 h (not shown) and extends into the eastern 
portion of the storm by 12 h (Fig. 8b). [In DRY144, the 
pattern is similar to DRY90, but the dry tongue has only 
just entered the storm at its southern edge by 12 h (not 
shown)]. The presence of dry air within the storm’s cir- 
culation produces two effects. First, it enhances the cold 
downdrafts in the southernmost portion of the precipi- 
tation immediately ahead of the dry tongue, which then 
expands to the south and east of the center. Second, it 


completely inhibits precipitating convection in the tongue 
of dry air. Differences in the radial distribution of azi- 
muthally averaged near-surface 6 e between DRY90 (Fig. 
9b) and CNTL are small and vary between positive and 
negative values through the first 3 days of simulation. The 
main exception is at radii greater than 90 km (i.e., the initial 
radial distance of the dry air) during the first 1.5 days of 
simulation, where DRY90 is generally about 1-2 K cooler. 

Very dramatic changes in near-surface 9 e occur in 
DRY0 (Fig. 8c) compared to CNTL and DRY90. Colder 
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Fig. 9. Radius vs time distributions of the azimuthally averaged 6 e (shading) at the lowest model level (90 m) for 
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interval is 2 K, starting at — 7 K, with differences <-5 K shown as thick contours. 


248 


JOURNAL OF THE ATMOSPHERIC SCIENCES 


Volume 69 


downdrafts lead to more expansive cold pools than in 
CNTL, particularly to the south and east during the first 
12 h. The surface pressure center is initially displaced 
toward the moist air during the first 12 h, but it then 
realigns with the circulation center and the core of 
higher 9 e air by 24 h (not shown). Figure 9c shows that 
the azimuthal mean 9 e pattern is much colder for DRYO 
than CNTL during the first 1.5 days, similar to CNTL 
between 1.5 and 2.5 days, and then lower than CNTL at 
later times. This latter result is partly because of 6 e in- 
creases in CNTL associated with its intensification, but it 
also results from cold downdrafts that are generated 
within DRYO (discussed below). The ~36-h length of 
time during which DRYO has colder 9 e values at the 
beginning of the simulation is comparable to the roughly 
2-day delay in intensification of the storm relative to 
CNTL. It is also evident from Figs. 5c and 8c that pre- 
cipitating convection is virtually always absent from the 
dry tongue during its initial intrusion into the pre- 
cipitation regions. This result suggests strong suppression 
of deep convection by the dry air through the combined 
effects of boundary layer cooling and entrainment of 
dry air into the rising cumulus towers. 

One might conclude from the above findings that the 
enhancement of cold downdrafts resulting from the dry 
intrusion is the primary cause for the delay in storm 
intensification in DRY144, DRY90, and DRYO, but the 
explanation appears to be more complex. For example, 
although DRY90 is weaker than CTRL, the differences 
in near-surface 9 e between CNTL and DRY90 are small 
(Figs. 9a,b) prior to day 4 when much of the delay in the 
intensification of DRY90 occurs. In addition, the rela- 
tively strong storm in DRYALL75 presents compelling 
evidence that cold downdrafts alone may not be the main 
mechanism for delay. This case produces the coldest 
downdrafts of the four simulations, with large areas 
at 12 h (Fig. 8d) with 9 e > 8 K cooler than in CNTL. 
The cold air spreads out symmetrically surrounding 
an inner region of moderate 9 e values (—344-345 K) 
that is similar to that found in the inner regions of all 
of the simulations. This simulation is colder than CNTL 
(Fig. 9d) at most radii outside the eyewall throughout the 
simulation and is up to 3 K cooler than DRYO at radii 
greater than 40 km at most times through the first 3 days 
of simulation. Despite the colder boundary layer air, 
DRYALL75 has an intensity that is nearly identical to 
CNTL throughout the simulation. 

Considering the inability of downdrafts to impede in- 
tensification in DRYALL75, combined with the strong 
wavenumber- 1 asymmetries in convection that are prom- 
inent in DRY144, DRY90, and DRYO, we suggest that the 
delayed intensification in these latter cases is strongly tied 
to their asymmetric structure. The asymmetric structure is 


a result both of the direct suppression of convection in 
the areas of intruding dry air (from lowering of boundary 
layer 9 e and entrainment) and of asymmetries possibly 
forced by asymmetric patterns of cold outflow that extend 
beyond the intruding dry air tongue. The remainder of 
this section will address this issue. 

To focus the discussion, Fig. 10 shows the time series 
of maximum azimuthally averaged tangential wind at the 
lowest model level (90 m) for the DRY270, DRY144, 
DRY90, and DRYO simulations. The intensity in the 
DRY270 simulation is similar to CNTL (see Fig. 1), with 
steady intensification from about 1 to 3 days, a period of 
more rapid intensification through day 4, and then mod- 
erate intensification through day 6. The DRY144 simu- 
lation begins similarly. By 2.3 days, it is about 2 m s~ x 
stronger than DRY270 but then undergoes a period of 
slow growth between 2.3 and 3 days. Between 3 and 
4 days, DRY144 undergoes moderate growth, but at 
about half the rate as in DRY270. As a result, the dif- 
ference between DRY270 and DRY144 becomes no- 
ticeable just after day 3 and increases through day 4, 
remaining steady thereafter up to day 6. In DRY90, the 
intensity also follows DRY270 up to about day 3. It then 
undergoes a 10-12-h period of weakening followed by 
slow strengthening through 4.4 days, rapid intensification 
through day 5, and then moderate intensification there- 
after. Finally, DRYO undergoes an approximately day- 
long delay in initial development, clearly tied to the 
general suppression of convection during the period of 
axisymmetrization of the dry air initially in the vortex 
core (see Fig. 7i). This delay is followed by a period of 
moderate intensification through about 3.5 days and 
then a period of no intensification through 4.5 days, 
after which moderate intensification resumes. Clearly, 
there are events that occur in the 2-4-day time frame 
that slow or prevent intensification in DRY144, DRY90, 
and DRYO. To further explain intensity differences, we 
next examine a characteristic of the convective evolution 
that is common to the three simulations and enhanced 
with increasing dry air. 

In general, the intensification of a vortex will be larger 
when the azimuthal mean diabatic heating occurs at a 
smaller radius closer to the high-vorticity core, whereas 
heating at larger radii away from the core will contribute 
much less to intensification (Hack and Schubert 1986; 
Montgomery and Enagonio 1998; Moller and Montgomery 
2000; Nolan et al. 2007a; Bui et al. 2009). In this regard, 
Fig. 11 shows the evolution of the radial distribution of 
azimuthally averaged diabatic heating at the 3-km level. 
As in Fig. 10, we use the DRY270 case as our reference 
for comparison because it has the same intensity as CNTL 
and includes the dry air. In DRY270 (Fig. 11a), the heating 
in the nascent eyewall increases gradually through the first 
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Fig. 10. Time series of maximum azimuthally averaged tangential wind speed for all simulations 
with an initial west-to-east-oriented dry air boundary for the first 6 days of simulation. Vertical 
lines indicate the start and end times for particular periods of slower growth compared to DRY270, 
with the line style matching the corresponding wind speed curve for the different experiments. 


3 days while contracting to smaller radius. This con- 
traction does not occur steadily, as in CNTL (not shown), 
but occurs with two brief periods of no contraction 
around 2-2.5 days and 3 days. The heating increases 
markedly with the onset of rapid intensification after 
3 days. A dashed line is used to indicate the axis of peak 
heating and is overlaid in the plots from other simula- 
tions to facilitate comparisons. Figure lib shows the 
pattern of diabatic heating for the DRY144 simulation, 
which is generally similar to DRY270 barring a couple 
of exceptions. First, prior to 2.3 days, the maximum 
heating contracts to smaller radius and is larger than in 
DRY270, consistent with the slightly greater maximum 
wind speed at that time. Just after 2.3 days, the heating 
peak stops contracting, the heating becomes smaller, 
and the rate of intensification slows. Second, a brief lull 
and radially outward displacement occurs near 3 days, 
followed by a more gradual increase in heating in the 
eyewall during the subsequent period of moderate in- 
tensification. The lull in and outward displacement of 
the heating exactly coincides with the separation of the 
tangential velocity curves in Fig. 10 at 3 days between 
DRY270 and DRY144. In the DRY90 simulation 
(Fig. 11c), the weakening phase coincides with a sharp 
outward displacement of the heating (indicated by the 
dotted curve) just after day 3, while the slow growth 
phase corresponds to a slow recovery of the radius of 
peak heating toward the DRY270 case. Finally, in the 
DRYO case (Fig. lid), the first period of no growth within 
the first 2 days corresponds to a general absence of 
significant heating as the dry air suppresses convection, 


while the second period of no growth beginning near 
3.5 days coincides with a large outward displacement of 
the heating relative to DRY270 between 3 and 4.5 days. 
Intensification resumes (see Fig. 10) once the heating shifts 
back toward smaller radius (near 30 km) at 4.5 days. 

The outward displacement of diabatic heating in 
DRY144, DRY90, and DRYO relative to DRY270 is 
related to the formation and outward movement of a 
highly asymmetric spiraling rainband that forms in each 
simulation. The asymmetry in precipitation increases 
when dry air begins nearer to the center of the vortex 
(e.g., section 4) and the subsequent rainbands encounter 
increasingly drier environments. For example, moving 
from DRY270 to DRY144, DRY90, and DRYO (Figs. 
12a-d), just radially outward from the rainband on the 
southwestern side of the storm at 3 days, 9 e is near 340, 
338, and 336 K, and then less than 336 K, respectively. 
As will be shown below, the interaction with their envi- 
ronment affects the duration and outward movement of 
the rainbands as well as the resulting cold-pool 6 e de- 
crease within the primary boundary layer inflow region. 

Figures 13-16 use 6-h averaged 2 near-surface 9 e and 
rain mixing ratio for the DRY270, DRY144, DRY90, 
and DRYO simulations to demonstrate the evolution of 


2 The fields have been smoothed to improve legibility of the 
figures. We use the Interactive Data Language’s (IDL) two- 
dimensional “median” filter of 7-point width. Median smooth- 
ing replaces each point with the median (the value with an equal 
number of values above and below it) of the two-dimensional 
neighborhood of a given width. 
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Fig. 11. Radius vs time distributions of the azimuthally averaged diabatic heating (shading) at 3 km for 
(a) DRY270, (b) DRY144, (c) DRY90, and (d) DRY0. The dashed line indicates the axis of peak heating in the 
DRY270 simulation. The dotted lines in (b)-(d) show the axis of peak heating in the other simulations where it lies at 
larger radius relative to the DRY270 case. The periods of slower growth (from Fig. 10) are indicated by the horizontal 
lines. 
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Fig. 12. Six-hour averaged S e (shading) at 3 km for the period 3-3.25 days from the 2-km grid for the (a) DRY270, 
(b) DRY144, (c) DRY90, and (d) DRYO runs. Contours show rain mixing ratios at 0.2 g kg~' interval (solid lines) 
plus an additional contour at 0.1 g kg -1 (dotted lines). The contour at 0.6 g kg -1 is highlighted by the thick contour. 


the rainbands that impact storm intensification in the 
2-4-day time frame. Dark red lines indicate the axes of 
peak inflow near the surface. In DRY270, at 2.5-2.75 days 
(Fig. 13a), the lowest 9 e air extends from the northern side 
around to the western and southern sides of the storm. 
Rain extends around the nascent eyewall, although it is 
weaker on the northeastern side. Over the next 12-36 h, 
the rain becomes increasingly asymmetric, with maximum 
values on the southern side and with the widest radial 
extent occurring between 3.5 and 4 days (Figs. 13c, d). 
During the same time, there is a gradual increase in 0 e 
values, particularly within about 100 km from the center. 
The degree of asymmetry in structure does not prevent 


rapid intensification during this time period since diabatic 
heating remains concentrated in the eyewall. 

In DRY 144, a prominent asymmetry in rainfall struc- 
ture is apparent at 2.5 days, with intense rainfall in the 
southern eyewall and little precipitation in the northern 
eyewall (Fig. 14a). At 3 days (Fig. 14b), precipitation 
increases on the southwestern side outside of the eye- 
wall, cooling the boundary layer along the main inflow 
region (red line). This time corresponds to the brief 
period of diminished and outwardly displaced average 
latent heating relative to DRY270 (Fig. lib). Over 
the next 24 h (Figs. 14c,d), as the rate of intensification 
increases, the rainfall in the eyewall intensifies and 
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Fig. 13. Six-hour averaged 6 e (shading) at the lowest model level (90 m) at the specified times from the 2-km grid 
for the DRY270 run. Contours show rain mixing ratios at 0.2 g kg -1 interval (solid lines) plus an additional contour 
at 0.1 g kg~' (dotted lines). The contour at 0.6 g kg -1 is highlighted by the thick contour. The red circle (of fixed size) 
is used as a reference to highlight changes in precipitation radius with time. Dark red lines are drawn along the axes of 
maximum radial velocities. 


contracts to smaller radius while 9 e values along and 
upstream of the peak inflow increase by ~3 K. 

The rainband development is even more dramatic in 
DRY90 from 3 to 3.5 days (Figs. 15b,c), where the highly 
asymmetric rainband spirals outward from the south- 
eastern eyewall around to the west. The rainband persists 
for —24-30 h, maintaining cooler 0 e values in the inflow 
layer and either weakening or significantly slowing in- 
tensification during this period. By —4.25 days (Fig. 15d), 
the spiraling rainband and cooler air has moved far out 
from the center while rain has reintensified at smaller 


radius in the eyewall and 9 e values have quickly in- 
creased behind the rainband. 

Finally, in DRYO, the very pronounced rainfall asym- 
metry evident at 2.75 days (Fig. 16a) gradually expands 
outward to larger radius on the southwestern side of 
the storm over the next 30-36 h (Figs. 16b-d), cooling 
the boundary layer inflow and substantially reducing the 
rain in the eyewall region. By 4.5-5 days and beyond 
(Figs. 16d-f), boundary layer 9 e values in the inner-core 
region behind the rainband increase quickly, intense 
rainfall redevelops in the eyewall, and the radius of the 
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Fig. 14. As in Fig. 13, but for the DRY144 simulation. 


eyewall contracts. Similar to other simulations, the ex- 
pansion of the outward spiraling rainband coincides with 
the period of no growth. 

In summary, as asymmetric dry air is initialized closer 
to the vortex core, it increasingly suppresses convec- 
tion, inducing an increasingly prominent wavenumber-1 
asymmetry in storm structure. After axisymmetrization, 
midlevel 6 e just outside of the moist core of the storm is 
lower when dry air is initially closer to the vortex center. 
With lower ambient 9 e values, the spiral rainband de- 
velopment becomes more pronounced and persistent, 
the 9 e in the area of strongest inflow becomes lower, there 
is a larger departure of latent heating from the vortex 
core, and there is a greater reduction in the intensification 


rate. While the relationships among rainband devel- 
opment, latent heating structure changes, and intensi- 
fication rates have been demonstrated, the exact means 
by which the dry air makes the rainbands more pro- 
nounced and prolonged and the extent of the role of 
boundary layer cooling remain unclear. 

6. Conclusions 

This study has focused on the impact of dry midlevel 
air on the development of tropical cyclones, specifically 
addressing the role of the dry air in enhancing cold 
downdraft activity and suppressing storm development. 
The WRF model is used to construct two sets of idealized 
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Fig. 15. As in Fig. 13, but for the DRY90 simulation. 


simulations of hurricane development with a non-SAL 
sounding: 1) a set of simulations with dry air located north 
of the vortex center by distances ranging from 0 to 
270 km and 2) a second set of simulations with dry air 
completely surrounding the vortex, but with moist en- 
velopes in the vortex core ranging in size from 0 to 
150 km in radius. 

For the first set of simulations, no impact of the dry air 
was seen for dry layers located more than 270 km north 
of the vortex center (~3 times the initial radius of 
maximum wind). As the dry air boundary was moved 
closer to the vortex center, the tangential flow of the 
vortex increasingly wrapped the dry air into the region 
of inner-core convection. The dry air suppressed initial 


convective development, leading to asymmetry of the 
convective vertical mass flux and slower storm de- 
velopment. Note that all simulations eventually reached 
the same steady-state intensity. 

For the second set of simulations, the presence of 
dry air throughout the domain, including the vortex 
center, substantially suppressed storm development, de- 
laying intensification at least 3 days. FLowever, dry air 
throughout the vortex is rather unrealistic. Observations 
suggest that most systems have a pocket or envelope of 
high humidity within the vortex core. When moist en- 
velopes (consisting of the non-SAL thermodynamic 
characteristics) were included within the vortex even out 
to a radius less than the initial radius of maximum wind. 




January 2012 


BRAUN ET AL. 


255 



a) 2.75-3 days 


-200 -100 0 100 
X (ton) 


200-200 -100 0 100 
X (ton) 


200-200 -100 0 100 
X (ton) 


200 


— I I I I 


345 348 351 354 357 360 363 K 


Fig. 16. As in Fig. 13, but for the DRYO simulation. 


the storm intensity evolved in a manner very similar to 
the control run without dry air, but the storm size was 
significantly reduced, consistent with the findings of Hill 
and Lackmann (2009). 

Downdraft cooling of the PBL during the initial in- 
gestion of dry air (the first 1-3 days) appears to have had 
little impact on storm intensification except in the 
DRYO simulation, where a clear suppression of de- 
velopment is seen during the first 2 days. After that 
time, differences in the rate of intensification are shown 
to be related to the development of an outwardly 
spiraling asymmetric rainband that increasingly slows 
storm growth as dry air is started closer to the storm 
center. Interestingly, DRYALL75 produces the most 
downdraft cooling and yet intensifies at a rate similar to 
CNTL, likely because the dry air maintains a symmetric 
zone of convection in the eyewall that locks the diabatic 
heating onto the region of higher vorticity in the storm 
core. This result suggests that strong downdraft cooling in 
and of itself does not necessarily inhibit intensification. 
Instead, it may need to be coupled with the production of 


storm asymmetries that have the potential to dramati- 
cally change the radial distribution of diabatic heating 
and push the heating peak away from the high vorticity 
core. 

The results above suggest that proximity of dry air 
near or even surrounding a moist vortex should not be 
interpreted as an indication of likely suppression of trop- 
ical storm development or a mechanism for storm weak- 
ening (e.g., Shu and Wu 2009). The dry air must approach 
very close to the inner core of the storm, to a distance 
comparable to or just outside of the radius of maximum 
winds, in order to slow the intensification of a developing 
tropical cyclone. Otherwise, the dry air apparently acts 
only to affect the size of the storm (Kimball 2006; Hill and 
Lackmann 2009). The results also suggest that the dry air 
must act in tandem with other processes such as vertical 
wind shear to produce more substantial inhibiting effects 
(e.g., Shelton and Molinari 2009). 

The simulations in this study involve highly idealized 
environments with no mean flow in the environment of 
the vortex. It is very possible that the impact of dry air 
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might be enhanced when combined with more complex 
environments, including steering flow and vertical wind 
shear. The idealized simulations obviously only repre- 
sent storms forming in conditions with little or no wind 
shear, but such conditions are not uncommon for storms 
that become major hurricanes. A future study will ad- 
dress the added complexity of mean flow and sheared 
environments. 
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